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5-H ' Abstract. We are concerned with global steady subsonic flows through general 

infinitely long nozzles for the full Euler equations. The problem is formulated 
as a boundary value problem in the unbounded domain for a nonlinear elliptic 
equation of second order in terms of the stream function. It is established that, 
when the oscillation of the entropy and Bernoulli functions at the upstream is 
sufficiently small in C 1 ' and the mass flux is in a suitable regime, there exists 
a unique global subsonic solution in a suitable class of general nozzles. The 
assumptions are required to prevent from the occurrence of supersonic bubbles 
inside the nozzles. The asymptotic behavior of subsonic flows at the downstream 
and upstream, as well as the critical mass flux, has been clarified. 



1. Introduction 

We are concerned with global steady subsonic flows through general infinitely 
long nozzles for the full Euler equations (without the isentropic and irrotational 
requirement). The two-dimensional steady full Euler equations take the following 
form: 

(f>u)xi + {PV) X2 = 0, (1.1) 

(pu 2 ) Xl + (puv) X2 +p Xl = 0, (1.2) 
(puv) Xl + {pv 2 ) X2 +p X2 = 0, (1.3) 

(pu(E + -)) X1 + (pv(E + -)) xa = 0, (1.4) 
P P 

where p, (u,v), p, and E denote the density, velocity, pressure, and total energy 
respectively. Moreover, 

!/-.2 , „.2n , P 



with adiabatic exponent 7 > 1. 
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Consider flows through an infinitely long nozzle given by 

= {(xi, x 2 ) ■ fi(x\) <x 2 < h{xi), -oo < X\ < oo}, 
with the nozzle wall d£l := W\ U W 2 , where 

Wi = {(xi,x 2 ) : x 2 = fifa), -oo < x\ < oo}, i = 1,2, 

as in Fig [U 



ft 



Figure 1. Infinite nozzle 



Suppose that Wi and W 2 satisfy 
/2(a?i) > 

/i(a?i) -> 0, f 2 (xx) 
fi(xi) ->■ a, / 2 (xi) 
and there exists a > such that 



II/, 



1 

• b> a 

< a 



for x\ G 
as xi — > 
as xi — )• 00 

* = 1,2, 



-00, 00 
00 in C 



in C 



2,a 
2,a 



(1.6) 
(1.7) 
(1.8) 

(1.9) 

for some positive constant C . It follows that £1 satisfies the uniform exterior sphere 
condition with some uniform radius r > 0. 

Suppose that the nozzle walls are solid so that the flow satisfies the slip boundary 
condition: 

(u,v)-n = ondn, (1.10) 

where n is the unit outward normal to the nozzle wall d£l. It follows from (II. ip and 
(fTTU]l that 



(pu, pv) ■ n dl = m 



11.11) 



for some constant m, where £ is any curve transversal to the x\— direction, and n is 
the normal of I in the positive x\— direction. 

If the flow is away from the vacuum state, it follows from fjl .21) — HI -41) that 

(u,v) ■ V(lnp- 7 hip) = 0, (1.12) 

which implies that ^ is a constant along each streamline, provided that the solution 
is C 1 -smooth. We assume that the entropy function is given in the upstream, i.e., 

IP 



(7 - l)pi 



S{x 2 ) 



as x\ 



-00, 



(1.13) 
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where S(x 2 ), defined on [0, 1], is the entropy. The sonic speed of the flow is defined 
by 




(1.14) 



By (JIT) and ([O]) . we obtain 



(u, «) • V(i(u 2 + v 2 ) + (^^) = 0. (1-15) 

This implies that \{u 2 + v 2 ) + which is called the Bernoulli function, is a 

constant along each streamline. We assume that the Bernoulli function is given in 
the upstream, i.e., 

(n 2 + v 2 ) ^p 

h t 7T > B(x 2 ) as xi -)• -oo, (1.16) 

where B(x 2 ) is defined on [0, 1]. 

Problem 1. Solve the full Euler system (jl.ip - (|1.4p wii/i i/ie boundary condition 
(jl.lOp . i/ie mass /Zux condition (jl.lip . and i/ie asymptotic conditions (|1.13p and 
(fTTT6|) . 



Set 



5= inf 5(x 2 ), 5= inf Bfe). 
z 2 e[o,i] x 2 g[o,i] 



The main results of this paper are the following. 

Theorem 1.1 (Main Theorem). Let the nozzle walls d£l satisfy (|1.6p - (jl.9p . and let 

S_ > and I? > 0. Then there exists 5q > suc/i i/iai, i/ 

\\(S-S,B-B)\\c^([0,l])<$ for0<5<5 , (1.17) 

and 

(S£^)'(0) > 0, (SB^)'(l) < 0, (1.18) 

there exists m > 25^^ such that, for any m € (5 1 / 4 ,m), there exists a global solution 
(i.e. a full Euler flow) (p,u,v,p) G C l,a {Cl) of Problem 1 such that 

(i) Subsonicity and positivity of the horizontal velocity: The flow is uniformly 
subsonic globally with positive horizontal velocity in the whole nozzle, i.e., 

sup(n 2 + v 2 - c 2 ) < 0, n>0 in U; (1.19) 

(ii) Far field behavior: The flow satisfies the following asymptotic behavior in 
the far fields: 

(a) As x\ — > —oo, 

j>-»Po>0, n-^n (x 2 )>0, (v,p) (0, po(x 2 ;p )), (1.20) 

Vp^O, Vn^ (0,u' (x 2 )), Vv^O, Vp -> (0, p' (x 2 ;p )) (1.21) 
uniformly for x 2 G i^i (<= (0, 1); 
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(b) As X\ — > OO, 

p^pi>0, u^ Ul {x 2 )>0, (v,p) (0,pi(x 2 ;pi)) ; (1.22) 
Vp^O, Vu^(0,u[(x 2 )), Vv^O, V P ^{0,p' x {x 2 - Pl )) (1.23) 

uniformly for x 2 £ K 2 (s (a,b), where po and pi are both positive constants, 

Po(X2;Po) = (( 7 -l)5(x 2 )) ' ^ l(x2;Pl)= (( 7 -l)5(x 2 )J ' 

and po, pi, uq(x 2 ) and u\(x 2 ) can be determined by m, S(x 2 ), B(x 2 ), and 
b — a uniquely; 

(iii) Uniqueness: The full Euler flow of Problem 1 satisfying (I1.19|) and the as- 
ymptotic behavior (11.20j) - (|1.23p is unique. 

(iv) Critical mass flux: m is the upper critical mass flux for the existence of 
subsonic flow in the following sense: either 

sup(u 2 + v 2 — c 2 ) — > as m — > rh, (1-24) 
n 

or there is no a > such that, for all m G (771,771 + a), there are full Euler 
flows of Problem 1 satisfying (|1.19p . the asymptotic behavior (|1.20p - (ll.23|) . 
and 

sup sup(c 2 - (u 2 + v 2 )) > 0. (1.25) 

mS(m,m+ir) Q 

The assumptions in Theorem 1.1 are required to prevent from the occurrence of 
supersonic bubbles inside the infinitely long nozzles. 

There has been some literature on the analysis of the infinitely nozzle problems. 
For potential flows, Chen-Feldman [5l [6] established the existence and stability 
of multidimensional transonic flows through an infinite nozzle of arbitrary cross- 
sections; also see Chen-Dafermos-Slemrod-Wang [S] and Kim p3]. Xie-Xin |16| 
established the existence of global subsonic isentropic flows and obtained the critical 
upper bound of mass flux under the assumption that the derivative of the Bernoulli 
function equals to zero on the two boundaries. For the steady full Euler equations, 
Chen-Chen-Feldmann [7] established the first existence of global transonic flows 
in two-dimensional infinite nozzles of slowly varying cross-sections; also see Chen 
[9]. Motivated by the earlier results, the focus of this paper is on the full Euler 
equations for the infinitely nozzle problem with general varying cross-sections by 
developing some useful new techniques. Some further related results can be found 
in Bae-Feldman [1], Canic-Keyfitz-Lieberman [4], Glimm-Ji-Li-Zhang-Zheng |13j . 
Serre [15], Yuan [17], and the references cited therein. 

We remark that the main difference between our results and those in |16] is that 
our results allow the varying entropy function, so that the far behavior of the density 
and the equation for the stream function is not only determined by the Bernoulli 
function, but also by the entropy function. Thus, it is not clear whether one can 
directly use the implicit function theorem to obtain the density with respect to the 
Bernoulli function at the upstream, which is the starting point of our study of this 
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problem. Furthermore, it is not direct to see how the maximum principle can be 
employed to locate these solutions of the stream function in the physical interval 
and then to extend the existence for small enough momenta which is obtained by 
the standard energy estimates to the critical mass flux only by simply making the 
assumption on the Bernoulli function on the boundary as in [16]. In this paper, for 
the far field behavior, we introduce the ratio of the entropy and Bernoulli function; 
then, by carefully defining the upper and lower bounds for the pressure, we find the 
far field behavior of the pressure with respect to the others. In order to use the 
maximum principle, we extend the entropy function via a special form, under which 
we find a condition on the ratio of the entropy and Bernoulli function with some 
power which looks like but not exactly the condition on the change of the momenta 
on the boundary at the upstream. With the uniform estimates from the maximum 
principle, we extend the existence of solutions to the critical mass flux. 

The organization of the paper is as follows. In Section 2, we reformulate the 
problem as Problem 2 by deriving the governing equation and boundary conditions 
for full Euler flows in terms of the stream function, provided that the Euler flow has 
a simple topological structure and satisfies the asymptotic behavior (|1.20p - (|1.23p . 
In Section 3, the existence of solutions to a modified elliptic problem is established. 
Subsequently, in Section 4, we analyze the asymptotic behavior of solutions in a 
larger class and show the uniqueness of the solution to the boundary value problem. 
This yields the existence of solutions of the boundary value problem for the stream 
functions. In Section 5, the existence and uniqueness of solutions of Problem 2 are 
established, and, in Section 6, some refined estimates for the stream function is 
derived. The proof of the main theorem (Theorem 1.1) except the existence part of 
the critical mass flux is provided based on the results of Sections 2-5. In Section 
8, we obtain the critical mass flux. Combining these estimates with the asymptotic 
behavior obtained in Section 4 yield the existence of full Euler flows which satisfy 
all the properties in Theorem 1.1. 

2. Reformulation of the Problem for Stream Functions 

In this section we introduce the stream functions for the two-dimensional steady 
compressible full Euler flows and derive an equivalent formulation for the full Euler 
flows in the nozzles. 

2.1. Equations. It follows from (jl.ip that there exists a stream function tp such 
that 

ipxi = ~pv, ipx 2 = pu- (2.1) 
Furthermore, from (I1.12p . we have 

P= 1 —S(^)p\ (2.2) 

7 

By (|1.15p . the Bernoulli law can be also written as 

l -\V^\ 2 + S^)p^ l =B^)p 2 . (2.3) 
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In the subsonic region, we have 

|V^| 2 <cV = (7-l)5(V)p 7+1 , 

which implies 

7 -i 25 

P 1 > 



(7 + 1)5' 

Let x = slW-l 2 and g(p,tp) = B(i>)p 2 - S(V>)p 7+1 - We obtain 



in the subsonic region. Hence, by the implicit function theorem, there exists a unique 
P = p{x-> '0) such that 

X = g( P , ^) = B^)p 2 - S{i,)p l+l . (2.4) 

From (|2.4p . we have 

1 _ B'p - S'pi 

Px ~ ~(~f + l)Sp-y -2Bp H ~ (7 + 1)5^-! -IB' 
Then we have 

PX! = PX^XX^XXXX + ^X2^X 2 X 1 ) + P^X! 

-i/'iiini - V^V^i + fpxAB' 'p 2 - 5' / o 7+1 ) 



( 7 + l)Sff - 2B P 

Px 2 = Pxi^ilpxixa + Tpx 2 1px 2 x 2 ) + Pipipx 2 

~^x^x^X2 - TpX2^X2X2 + 4>x 2 (B'p 2 - 5'p 7+1 ) 



(2.5) 



(2.6) 



(7 + l)Sff - 2Bp 

Now we reduce the Euler system into a second-order nonlinear equation. Multi- 
plying equation (|1.3p by (7 + l)5p 7 — 2Bp, using expressions (|2.5p - (|2.6p . and making 
algebraic manipulations, we obtain 

K (ay ty, Vi>)^ XiXj - Fty, V^)) = 0, (2.7) 

where 

on(^,VV) = (7- l)«V +1 -<, 
ai 2 (^, V-0) = a 2 i(ip, V-0) = ip Xl tpx 2 , 
a 22 (^V^) = ( 7 - 

W) = ^^ / o 7+3 (7cSi3' - 25'B + cScS'p 7 " 1 ) • 

7 

If u > in f2, then t/> X2 > 0. Thus we have 

^,V^ XiXj =F(^,V^). (2.8) 
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Multiplying by ((7 + l)5/> 7+3 - 2Bp 4 ) , we obtain 

V-( — ) = B'p--S'p 1 . (2.9) 
P 7 

In summary, we have the following proposition. 

Proposition 2.1. For any smooth flow away from the vacuum state in the nozzle 
Q satisfying (|1.7p - (|1.9p . if the flow is globally subsonic and 

u>0 inU. (2.10) 

Then the new system formed by (|2.1|) - (|2.3p and (|2.9|) is equivalent to the original 
Euler equations (ll.ip ~ p.4p . 

The previous derivation is obviously invertible for the subsonic flow, so we omit 
the details of the proof for Proposition 12.11 In order to establish the existence of 
solutions to system (|l.ip - (|1.4p . it suffices to establish the existence of solutions to 
system (|23 ]) -(|23 ]1 and CUD satisfying (^TO]) . 



2.2. Relations between S(tp), B(tp), and the Asymptotic Behavior of ip at 

xi — > ±00. First, it follows from (jl.lOp that the nozzle walls are streamlines, so ip 
is constant on each wall. By (jl.lip and the fact that tp X2 > since u > 0, we have 

< tjj < m in 0, ip = on W\, ip = m on Wi. (2-H) 



Then we study the density-speed relation by using the entropy relation (12.21) and 
the Bernoulli law (|2.3p . Here, unlike isentropic flow which does not need to study 
the entropy relation, i.e., the entropy function S(x2), we start from the ratio of these 
two functions B and S as follows. 

Let D(x2) = (BS t)(x2). For any s > 0, p(s) = -L- s1 -i > is the unique 
solution of 



7P(g)^ 
V 



7 



Moreover, from (|2.3p . the speed 



q(p,x 2 ;s) = y2^(x 2 )( S -(-^ T ) 7 "). 



Hence, for fixed s and x 2 G [0, 1], q is a strictly decreasing function of p on [0, p(s)]. 
By the definition of p(s), one has 

q(p(s),x 2 ;s) = < c(p(s),x 2 ). 

Now we claim that q(0,x 2 ] s) > c(0,x 2 ). Indeed, 



q(p, x 2 ; s) — > Y 25*t (^2)s > as p — > 0, 

and, by the definition of sonic speed, c(0,x 2 ) = 0. Thus, q(0,x 2 ) > = c(0, x 2 ) 
This completes the claim. 
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l 7—1 

Since c 2 (p,x 2 ) = (7 — l)S{x2)p' 1 ~ 1 = (7 7_1 (7 — l)S{x 2 )) 1 p~ is an increasing 
function of p, there exists a unique p(s) G [0, p(s)] such that 

c 2 (p(s),£ 2 ) = q 2 (p(s),x 2 ;s). 

More precisely, 

In summary, we have 
Lemma 2.2. There exist p = p(s), p = p(s), a^e? r = T(s,X2) suc/i £/iai 

^(^(fiMj^^W, (2 . 12) 
s ^, )( ^)^ + I^p> =B(l2) , (2 . 13 , 

c 2 (p(s),i 2 ) =r 2 (s,i 2 ), (2.14) 

where p(s), p(s), and r(s,X2) are i/ie maximum pressure, critical pressure, and 
critical speed, respectively, for the fixed ratio s of the Bernoulli function and the 
entropy function. 

Then direct calculations show that 

#>0, ^>0. 
as as 

Clearly, p(s) < p(s) for s > 0. By the continuity and monotonicity of p(s) and p(s), 
there exists a unique 5 > such that 

p(D + S)=p(D)- (2-15) 

where D = inf D(x 2 ). 

~ 12 e [o,i] 

Moreover, it follows from ()2. 15j) that there exists a uniform constant C > such 

that 

' C" 1 < p(D) < p(D) = p(D + 5) < C, 

C- 1 < p'(s) < C, C7" 1 < p'(s) < C itse(D,D + S), ( 2 .16) 

C7" 1 < S^(x 2 ) (^k) 2 ^ <C if p G (j>(2), *>(£ + £))■ 

Hereafter, C denotes a generic constant which depends only essentially on S" and B. 
If 5(x2) and B{x 2 ) satisfy 

\\(S-S,B-B)\\c^([o,i])<S, (2.17) 

then 

5= sup J^A- <D + C5. (2.18) 

«26[0,1] 5 ,i /7( 2 ; 2 ) 

Finally, we study the behavior of S and B at the upstream and downstream in the 
far fields of the nozzle where the flow may have certain simple structure. Indeed, for 
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the flows satisfying ()1.20|) (|1.25|) . one can determine po, pi, po(x2), p±(x2), uo(x2), 
and 1/1(3:2) first. 

If the flow satisfies (11.20p . then 



+ S{x 2 )pl 1 {x 2 ;po) = B(x 2 ), u (x 2 ) > 0, po(x 2 ;po) = ( 



2 V(7-l)S(x 2 ), 

(2.19) 

and 

/ Po{x 2 ;p )u (x 2 )dx 2 = m, (2.20) 
Jo 

which imply that 

u (x 2 ) = ^2(B(x2)-S(x2)pl- 1 (x 2 ;p )), (2.21) 

and 

m = j p (x 2 ;po)\/2(B(x2) - S(x2)po~ 1 (x2;po))dx 2 . (2.22) 

Lemma 2.3. Let 5 < §. It follows from (pTL5]l that p(D(x 2 )) < p(D) < p(D). Then 
we have 

(i) For given S(x 2 ), B(x 2 ), and m > ; (|2.22p has a solution po satisfying 
po € (p(D),p(D_)) such that po(x2',po) and uo(x2) > satisfy (|2.19p for 
x 2 €[0,l]; 

(h) If S > and \\B — jB||c«i,i ([o,i]) = ^ — f or some small 5q, then there is a 
positive constant C such that 

'C- X W <p(D)-Po<C, 
C-HV < u (x 2 ) < C, (2.23) 

The proof of this lemma is as follows. 

Result (i) is for obtaining a global subsonic flow in the nozzle. Clearly, from (|2.19p 
-^{f Po(x 2 ; P o))^2{B(x2) - S(x 2 )p2- 1 (x 2 ;p ))dx 2 ) < 0, ^|^) > 
forpo G (p(D),p(R)). Let 
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It follows from ([27161) and (f2"T7D that 

£ 0(D^2)^2S 1 h(x 2 )(D(x 2 ) - (S 1 hg(D;x 2 )p~ 1 )dx 2 

= J* Q(D\x 2 )yf2S 1 /'r(x 2 )(D(x2) - D)dx 2 
< Cb 1 ' 2 . 

In addition, 

J Q(D; x 2 )^2S l h(x 2 )(D(x 2 ) - (S l h (x 2 )g(D; x 2 ))^)dx 2 

> J g(D; x 2 )^2S^h{x 2 ){D- {S^h e {D-x 2 ))^)dx 2 
= J q(D; x 2 )^2S(x 2 )(e(D; x 2 ^~ l - e (D;x 2 p~ 1 )dx 2 

= J g(D; x 2 )^2S(x 2 )(g(D + 5; x 2 )^ - q(D; x 2 )^)dx 2 

> J q(D; x 2 )^2S(x 2 )(e(D + 5;x 2 p- 1 - g(D + 5/2; x 2 )^ l )dx 2 

Therefore, for any (3 G (0,1/3), there exists So G (0,5/2) such that f|2.22|) has a 
unique solution po G (p(D),p(D_)), if < 5 < So and m G (5^,mi), where mi 
satisfies C" 1 ^ 1 / 2 ; m-i > 25q > C6 1 / 2 . Later on, for simplicity, we will choose 
/3 = 1/4. However, all the results hold for /3 G (0, 1/3). 
By virtue of (|2.22p . one has 

m = J Q l P o{x 2 ;po)^2S l h{x 2 )(D(x 2 ) - (S 1 h(x 2 )po(x 2 ;p )p- 1 )dx 2 

= ! Q X Po{x 2 ;po)^2S 1 h(x 2 ){D{x 2 ) -D + D - {S^h (x 2 ) P o(x 2 ;po)V~ l )dx 2 

< CJ 1 P o(x 2 ;poWS+(S 1 ^(x 2 )q(D;x2)P- 1 - {S 1 h{x 2 ) P o{x 2 ;po)V- 1 dx 2 . 
Thus, we have 

S + (S^ix^giD;^))^ 1 - {S lh (x 2 )po{x 2 ;po)V~ l > C'H 2 ? . 
Since (3 < 1/3, then there exists So G (0, 5o) such that, if < S < So, then 

m^y- 1 - p y - 1 (x 2 ;p ) > c-h 2 ?. 

Consequently, if S_ > and — ^S||c?i,i([o,il) = ^ < there exists a constant C > 
such that (]2.23p holds. This completes the proof of Lemma 12.31 

Next, to determine the asymptotic states in the downstream, we parameterize 
the streamlines in the downstream by their positions in the upstream. Using (jl.20p . 
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(fL22l . and (l27TOl) . we define 

y = y{s) for AG [0,1], (2.24) 

such that 

S{s)pl\s-p Q ) + ^ = S(s) P r 1 (s; Po ) + Ul (y(s)) > 0, (2.25) 

J* p (t;po)u (t)dt = f* is) Pl (t; Pl ) Ul (t)dt, Pl (s; Pi ) := ( (t J^ (s) ) 7 , (2.26) 

2/(0) = a, y(l) = 6. (2.27) 

Then the streamline which starts at (— oo,s) ends at (oo,y(s)). 

The next procedure is similar as before, where we consider pi instead of pi for the 
monotone relationship between them and for simplicity by recalling 
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6{D]X2) := V (7-D^2) J ^ g " (AX2):= l ( 7 -l)5(x 2 ) J • 

The mapping in (I2.24p is well-defined due to condition (|2.25j) and (I2.26j) . In fact, 
({2T26D deduces that 

p {s;p )u (s) = pi(s;p 1 )ui(y(s))y'(s). (2.28) 
If pi < g(D.;x2) < e(D;x2), then 

5( S )(^- 1 ( S ;p )-pr 1 ( S ;Pi))+^ = B(s)-S(s) P r\s) > B^-S^^s) = . 
Then we have 

' ffy _ po(s;po)u (s) 

ds 



pi(s;pi)^/2S(s)02 1 {s;p a )-p1 1 (s;pi))+u%(s) ' (2.29) 
.2/(0) = a, 

where the pressure in the downstream pi satisfies 

fl>(«;K>)"o(«) & = 6 - a. (2.30) 

^ Pi( S ;pi)V2S( S )(pr 1 ( S ;Po)-pr 1 (^;Pi)) + n2( s ) 

It remains to show that there exists p\ G (p(D), p(D_)) satisfying (j2.30p . As in the 
proof to Lemma [231 we find that, for p\ G (p(l)), p(L>)), 

d f 1 Po(.s;p )u {s) 



dpi Jo 



Pi(s;pi)^2S(s)(p2 1 {s;p )- P J 1 (s;p 1 )) + u§(s) 
On one hand, there exists 5o G (0, So) such that, if 5 < So, then 

r 1 p (s;p )u (s) ^ s 

m\s)\J*S{s)(pZ-\s-,p )-e{&,s)^)+ul{s) 
_ f 1 po(s;po)u (s) ^ 

g(D;s)^/2S ,1 /7( s )(D(s)- ; D) 



ds > 0. 



> c^- 1 )/ 2 > b 



a. 
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On the other hand, 

r 1 po(s;po)u (s) ^ s 

_ r 1 po{s;p )u (s) ^ s 



< 



e (D;s)V25 1 /7( s )( jD ( S )-(5 1 /7p(Z)))7-i) 
I 1 
7^1 ) t rl (s)p (s;p )u (s) j s 

-yp(D)) JO \/2(D-{S 1 /~ip(D))~<- 1 ) 

( 7-1 N \ 7 rl 5^7 (s)po(s;p )n (s) 



fe^ /0 ;2((^l)^-(^))^)^ (2-31) 
< ( " m _ max^efo n (S%(x 2 j) 

- y ~ IPiD) \f(^)^-c-^)^) a[,u ] 

<gf<6-a. 7 

Thus, there exists a unique p\ G (p(.D),p(.D)) such that (j2.30p holds, provided 
that < 5 < 5o and m G (<5^, 777,2) for some <5o small enough and 26% < m 2 < 
min{mi,C _1 (6 — a)^ 1 / 2 }. Once pi determined, y(s), pi(s;pi), and ui(s) can be 
obtained from (|2.25j) . (|2.26p . and (|2.28p . Therefore, the above calculations yield the 
following proposition. 

Proposition 2.4. Let S, B > 0, and D be the ratio of S 1 ^ and B. There exists 
5q > such that, for any S, B G C 1,1 ([0, 1]) satisfying (|2.17|) with 5 <5o respectively, 
there exists rh > 25q , (3 G (0, |), such that 

(i) Existence: There exists solutions (po,Uo,Po) to (|2. 19j) f)2.20|) and (pi,Ui,pi) 
to ^M-^M if™ e (S^fh) mft pJiyfPj) = ^r^-yj = 0,1, y = x 2 
and yi = y(s); 

(ii) Subsonicity: po, p\ G (p(D),p(D)) ; 

(iii) Limiting behaviors: Either po —> p(D) or pi — > p(D) as m — >■ rh, 

where D = sup D(x 2 ) and 
x 2 e[o,i] 

m = sup{s : m G (S s) such that there exist po, p\ G (p(D),p(D))}. (2.32) 

Proof. Results (i)-(ii) are direct corollaries of Lemmas I2.2H2.31 It suffices to verify 
(iii). 

For m G (8",m 2 ), Po,Pi G (p(D),p(D_)). For fixed S(x 2 ) and B(x 2 ), po decreases 
as m increases. If 



m — > m 

then 



j\(D;x 2 )^2S l h(x 2 )(D(x 2 ) - (S l hg(D))^)dx 2 , 
Po->p(D). 
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For fh defined in (|2.32p . fh G [m2, fh). Note that both po an d pi are uniformly away 
from p(D)- If neither po nor p\ approaches to p(D) as m — > fh, then there always 
exist po,Pi G (p(D),p(D_)) for m G (5^,fh + e) for some small positive e, which 
contradicts with the definition of fh. This completes the proof. □ 

2.3. Reformulation of Problem 1: Problem 2. Let X2 be the coordinate in the 
upstream. Since po(^2;Po)^o(^2) > for X2 G [0, 1], ip is an increasing function of 
X2. Thus, we can represent X2 as a function of ip, which is defined by 

X 2 = k(V>), < V < m, 

so that 

ip{X 2 ) = / p (s;p )u (s)ds. 



Jo 

It follows from Proposition 2.1 that, if (|2. 10[) holds in SI, through each point 
(xijX?) G fi, there exists a unique streamline which starts from the upstream. 
Along each streamline, the stream function is a constant by the definition. Therefore, 
through any (x%, X2) in the nozzle, there exists a unique streamline from (—00, 
with tjj = ip{x\,X2). Thus, we denote 

S = S(k(i/>)), B = B(K(tp)) for < ijj < m. 

Then our main task in the rest of the paper is to solve the following problem: 

Problem 2 (Reformulation of Problem 1). Seek a solution of the boundary 
value problem: 

i v • (^mh)) = WMIWf , V) - iS'(V^(|VVf in SI, 

\y ~ h{xi)-h(xi) m on 
such that 

(i) The flow field induced by 

p = p{\VM\i>), u = ^, v = -^, p =l^±s(^ 

P P 7 

satisfies (TOUD - (fL25l) ; 

(ii) In the upstream, we have 

rX 2 

xl>{X 2 )= / p (s;p )u (s)ds, 0<^<m. (2.34) 
J o 

3. Existence of Solutions of a Modified Boundary Value Problem 

There are three main difficulties to solve the boundary value problem (|2.33|) . The 
first is that equation (|2.33p may degenerate at the sonic states. The second is that, 
although the entropy and Bernoulli function are well-defined on [0, m], the density 
p is not well-defined for arbitrary ip and IV^I- The last is that the problem is in 
an unbounded domain. Our basic strategy is to extend the definition of S{ip) and 
B(tjj) appropriately, introduce the elliptic cut-off to truncate |V^| in p(\X7ip\ 2 , ip) in 



11 
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a suitable way, and use a sequence of bounded domains and solve the problems on 
it to approximate the original one. 

In this section we first introduce a modified problem and then solve it, which can 
be indeed used to solve the original problem with the asymptotic behavior in §4. 

Set 



a(s) 



and 



b(s) 



fs'00 

S'(0)*±p 
10, 



(&)W 



We define 



if < s < m, 

if m < s < 2m, 

if — m < s < 0, 

if s > 2m or s < — m, 

if < s < m, 

if m < s < 2m, 

if — m < s < 0, 

if ip > 2m, or s < —m. 



f B(0) 



S(s) = S(0) + ^ a(t)dt, B(s) = &(s) + J Q b{t)dt) . (3.1) 

Then (5, B) G C 1 ' 1 (IR). We remark here that the definition of b(s) in this particular 
form instead of B' itself is for some technical reason, roughly speaking, due to the 
maximum principle. Moreover, since m > 8@, there exists a suitably small <5i such 
that, when 8 <6\, 

<B-C5 <\ul(s) + S{s)pl~ l < sup B(x 2 ) + C5, u {s) > 

for some C > 0, where ||(5 - 5(0), B - JB(0))|| c i,i( K ) < 

In the rest of the paper, we will always use the following notations: 



It is easy to see that 



9|V^| 2 



Pi(|V^| 2 ,^) 



ap(|v^l 2 ,^) 

dip 



1 



2p(c 2 - jV^IVp 2 ) 

goes to — oo when the flow approaches the sonic state from the subsonic states. To 
avoid it, we introduce the following cut-off function. For e > 0, let 



CoOO 



if s < -2e, 



(3.2) 
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be a smooth increasing function such that |£q| < 1- We define 

A 2 (|VVf , V) := CodWf - (7 - 1 W)P 7+1 ) + (7 - 1W)P 7+1 , (3.3) 

where 

iA 2 (|V^| 2 ^)+5(V0/^ +1 =m)f- (3-4) 
A direct calculation shows 

Sij{q,z) = p(\q\ 2 ,z)6ij - 2p 1 (\qf,z)£ i £j, (3.5) 
and ^ 

pi{\V4>\ 2 ,^) = ^ = < 0. 

Obviously, there exist two positive constants A(e) and A(e) such that 

A|£| 2 <,%((/, z)£^-<A|£| 2 (3.6) 

for any z G M, g € M 2 , and £ € M 2 , which means that the modified equation is 
uniformly elliptic. Thus, instead of solving Problem 2, we first solve the following 
problem: 

Problem 3 (Modified Problem). Seek a solution to the boundary value problem: 

-f( ) ( 3 - 7 ) 
such that 1 1 ip has a uniform upper bound. 

Proposition 3.1. Let the boundary d£l satisfy (jl.6p - (jl.9p . Then there exists < 
Si < min{5o,5i} such that, if \\(S — S_, B — 5)||c7i.i([o,i]) ^ ^ < 5 < <5i and 
m € (6^, mi) with mi = 25 /3 / 2 < m, where m is defined in Proposition 2.2, then 
Problem 3 has a solution ip G C 2,a (Q) satisfying 

M<C(e,J), |V^| 2 <( 7 - W +1 -2e (3.8) 
for some e > 0, so | Wf < S(e) - 2e with^e) := (7 + 1) (ff + J + e) ( ff^sj) ) ^ ■ 

Proof. The proof of the existence part is standard via approximation by the corre- 
sponding problems on bounded domains, while inequality (|3.8I) is crucial here, since 
S(e) depends not only on B but also on S 1 for the non-isentropic flows. We divide 
the proof into four steps. 

1. First, we use a sequence of boundary value problems on bounded domains to 
approximate Problem 3 on the unbounded domain. Since the key point is to obtain 
estimate (|3.8p . we focus on the following boundary value problem: 

V . { m = B'p-lS'p inn L , 

t( , (3-9) 
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where Ql satisfies 

{(xi,x 2 ) : (xi,x 2 ) G £1, \xi\ < L} C £Il C {(xi,x 2 ) ■ (x\,x 2 ) & fi, \xi\ < 4L} 

for all positive constants L > Lq > 0, with Lq sufficiently large, and G C 2 '" 1 , < 
ol\ < a, satisfies the uniform exterior sphere condition with uniform radius ro, 
< r < r. 

2. Equation (j3.7|) can be written as 

" P 2 |W>| 2 = (B' - -S'p*- 1 )?, (3.10) 

7 

where the repeated index is the summation with respect to the index from now on 
and 

- _ -2^p-( 7 -l)C^ + (7 + l)^ f3n) 
' AB - ( 7 + lYSp^ + ( 7 2 - V" 1 ' 
Therefore, f|3. lOj) becomes 

SijiV^^d^ = .F(V^), (3.12) 

where 

J-(V^, VO = {& - -S'p- l )p z + ,5 2 |VVf . 

7 

Instead of (|3.9p . we first solve the following problem: 

' A ij (V'il>,i>)d ij i> = T(yi>^) inO L , 



/ _ s2-/iQn) an 
w ~ f2{x 1 )-f 1 (x 1 ) m on cni L , 



(3.13) 



where -F(V - 0,^) = (-B' — -5' 1 )p 3 + P2A 2 for dealing with the fact that .F has 
quadratic growth in |Vi/>|. By the definition of £, 5, and ,6, we have 

|^(V^,^)| < Ctf. (3.14) 

3. Then, by the standard existence theory of elliptic equations, there exists a 
solution tpL to (|3.13p . Furthermore, writing ipj^ = m.m{tpL, 0} and ip^ = niaxj^, 0}, 
by the maximum principle with the source term (cf. Theorem 3.7 in 

C C 
min — - sup |.F| < ipi < sup + — sup |.F| , (3.15) 

A q l qq l A q l 



where C = e d — 1 with d = sup{f 2 (xi) — f\(x\)}. Then we have 

_ CS i-W _ C gi-P < t f Jk < m + C5 l ~ 2 P + OS 1 ' 13 for k sufficiently large. 

Moreover, one can obtain some other estimates for ipk. In fact, we can use the 
following more precise form with the same notations and symbols as those in Chapter 
12 in H2, 



[u]i, a < C( 7 , n) (l + \\Vu\\ + ^) . (3.16) 
Here, C( 7 , fl) depends only on diam(O) and the C 2 -norm of dQ. 
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Applying estimate (|3.16p to problem (|3.13p deduces that there exists fx = > 
such that, for any xo £ and for ip^ with k > 4L, we have 

&k}i^ Bl ( X0 )nn L < C(j, \\hh, ||/ 2 || 2 )(l + ||V^||o ;Bl (x )nc L + (3-17) 

Furthermore, using the interpolation inequality and the maximum principle (|3.15p . 
we obtain 

II^Jfe||l;Bi(*o)nni < ll/ilb. \\hh){lHy^k\\o-,B 1 (x )nn L + ^^-)+C v {m+^^-), 

where C > is the same constants as that in (13. 15H . Taking rjo sufficiently small so 
that nC(§, \\fxh, II/2II2) < \ if V < Vo, then 

\m\i;B l{xo )nn L < vC(j, H/1H2, ||/ 2 || 2 )(1 + + C„(m + (3.18) 

Thus, the Holder estimate (|3.17j) becomes 
IIV'felli^Bi (lo)nsJL 

= ||V'fc||l;Bi(xo)nn£ + [^fc]l,At;Bi(x )nn i 

< (1 + C7(A ll/!^, ||/ 2 || 2 ))||^|| 1 .B 1 (x )nn £ + C(§, \\hh, ||/ 2 || 2 )(m + 
<C7(A || /l || 2 ,||/ 2 || 2 )(i + m+ &). 

(3.19) 

Since, for any x,y £ Ql, 

\Vip k (x) - V-0fe(y)| ( IIV'fc||i, At ;Bi(a : o)nn i if y e -81(2:0) n n L , 

\ x ~ f I" " l2||^||i;B l{ xo)n^ if y £ Bi(a;o) H «l, 
which, together with f|3. 18|) and (|3.19p . yields the following Holder estimate: 

[AW-,n L < C(j, H/ilb, ||/ 2 || 2 )(1 + m+ (3.20) 
Thus, it follows from the standard Schauder estimate that 



lj kh,a;B 1/2 (x )nn L < C(p ll/illc* 2 ^) ll/2||c 2 ."i m 5 - ^° )- 



Thus, 



H^fclka;^ < <?(-, ||/i|| 2 , a , ||/ 2 || 2 ,a, ""i, ^T"^)- ( 3 -21) 

4. Using the Arzela-Ascoli lemma and a diagonal procedure, we see that there 
exists a subsequence ^ such that 

V'fc; — )• ^ i n C 2, ^(K) for any compact set if C and ?? < a. 
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Here, ip satisfies the following problem: 

'iy(V^,VMi^ = ^(V^,V) inn, 

X2-fl{x 1 ] 

/a(xij-/i(ii)' 

with the estimate 



1> = f.rffii "' on ^n, 



L;fi < 7/C(A, ||/i|| 2 , ||/ 2 || 2 )(1 + + C^TB + ^) 

< 7?C(A, HMI2, ||/ 2 || 2 )(1 + C8) +C v {m + CS), 

where r\ G (0,770) and C depends only on #o, m, A, and A. Next, we prove that 

|V?/f < (7 - l)5p 7+1 - 2e. 

Otherwise, 

( 7 -l)5pT +1 < |Vyf + 2e 

< 7 ? C(A,||/ i || 2 )(l + C5) + C^(?n + C7 ( 5) 



(3.22) 



< (7"1)5(- 



25 > ^ 



7-1 



'(7 + 1)5 
Thus, 

/ 2g x^i 

and 

2V - (7 + l)Sp< +1 > 0, 
which contradict with the fact that 

2Bp 2 - (7 + W +1 = Co(|V^| 2 - (7 - < 0. 

Thus, the solution ij) satisfies 

IV-^I 2 < S(e) - 2e (3.23) 

for any 5 G (0,«5i) and m G (^ , 2jf /2 ) . Then follows from (^22]) and ([333]) . 
Furthermore, (|3.20p and ()3.21[) yield the following higher order estimates 

ll^lli, M; a<C(pll/ill2,||/2|| 2 )(l + m+&), (3.24) 

and 

U\kn<C(j,\\fih,aA\f2h,a,m,^). (3.25) 

This completes the proof. □ 

Remark 3.2. Estimate (|3.8p in Proposition 13.11 implies that the cut-off function 
introduced in (13.21) and (13.31) can be removed. 



GLOBAL STEADY SUBSONIC EULER FLOWS THROUGH INFINITELY LONG NOZZLES 19 



4. Far Field Behavior of Solutions of Problem 3 



In this section, we study the far field behavior of solutions to Problem 3. We now 
show that the solutions to Problem 3 satisfy the asymptotic behavior (|1.20p - (|1.25p . 
and < ip < m. From this, we can remove both the extension and the elliptic cut-off 
(|3.7p . Therefore, these solutions solve Problem 2. In addition, the stream function 
formulation is consistent with the formulation of Problem 1 for the non-isentropic 
Euler system in the infinitely long nozzle, as long as the flow induced by a solution 
to Problem 2 satisfies (|1.20p - (jl.25p and ()2.10p . Furthermore, the far field behavior 
is crucial also for the consequent result of the uniqueness of the solutions. First we 
have 

Lemma 4.1. For e > 0, there exists 5 2 £ (0, Sq] such that, if 

(i) 11(5-5,5-5)1101,1 <5<5 2 ; 

(ii) to G (8",fh), where fh is defined as Proposition [27 



then there exists a function tp that satisfies 

tp — > tp as xi 

and 

X2 



-00, 



tp(xi,x 2 ) = tp(x 2 ) = / p {s;p )u (s)ds, 
Jo 



(4.1) 



where po and uq are uniquely determined by S, B, and m as §5, so tp is independent 
of Xl. 

Proof. The proof is based on the blowup argument in combination with the energy 
estimate, which consists of three parts: The first is for the existence of tp, the second 
is for the independence of ip of x±, and the third is the explicit form (|4.ip for ip. 

1 . Existence of the far field function. It is convenient to introduce a new coordi- 
nate to flatten the boundary walls of the nozzle, as follows: 

(ti(x 1 ,x 2 ) = Xl, 
\t 2 (xi,x 2 )- X2 ~ h ^ 



f2(xi)-fl(xi) ' 



then the nozzle becomes (—00,00) x [0,1]. Obviously, the coordinate transform is 
reversible, since 



det 



- dti_ 


dti - 






dx\ 


8x2 






9*2 


dti 






. 8X! 


8x2 - 











1 



' /2(zi)-/l(zi) 



f2(xi) ~ fl(xi] 



^0. 
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In addition, we remark here that the equation does not change the type of ellipticity 
under the coordinate transformation, since 

a . _/aii(x 2 -/i(si))(/2(a?i) ai2 \ a , 

= a u dt ltl ip - 21 — — — —= j— — — — -Y )d tlt2 ip 

V (/2(xi) - /i(xi)) 2 (/ 2 (xi)-/i(x!))/ 

| ^ii(*2 - /i(^i)) 2 )(/2(^i) " /i(^i)) 2 2a 12 (x 2 - AOn)) (/^(xO - /{On)) 



(/2(xi)-/i(x!)) 4 (/ 2 (zi)-/i(*i)) 3 



«22 



(^(Xl)-/!^!)) 2 



+ lower terms (involving d^ip and 

In the new coordinates, define 

?// n ) = ip(ti(xi - n,x 2 ),h(xi - n,x 2 )). 

For any compact set K C (—00,00) x [0,1], it follows from (|3.25p and the C 2,a - 
bounds of the walls f\ and f 2 that 

Hc^rR:) < C for n sufficiently large. 

Then, as in Step 4 of the proof to Proposition 13.11 there exists a subsequence 
such that 

^(n k ) ^ in ^(fT) (4.2) 

for any compact set .FT C (—00, 00) x [0, 1] and any $ € (0, a). From (jl.6p - (|1.9p and 
(|3.25p . and the facts that fi(x\) — > and f 2 {x\) — > 1 in C 2,a as xi — >■ —00, which 
also means that //(sgi) — > in C 1 '" as x\ — > —00, then satisfies 

V • (^Jrjj) = (S'p- i5'^)(|V^| 2 ,^) in I?, 
^ = on x 2 = 0, (4.3) 
■ip = m on x 2 = 1, 
where D = (—00,00) x (0, 1), and tp also satisfies 

\i>\<C(e,6), |V^| 2 <E(e)-2e. (4.4) 
Thus, by similar arguments as in §3, on any compact set E C (—00,00) x [0, 1], 

\m C ^(E) <C(e,S), 

and 

U\\ C 2, a{E) <C(e,5). (4.5) 
Thus, ip G C 2 ' a (D). This completes the first part. 
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2. Differentiate the equation in (|4.3p with respect to x\ and set u = ip X i- Then 

*(^>")-*(^^)=^«»^*)«*'. 

( 4 - 6 ) 

where Aij(q,z), @(q,z), and $(q,z) satisfy 

Aij(q, z) = p{\q\ 2 , z)8ij - 2pi(\q\ 2 , z)qtqj, 

Q(s,z) = B"(z)p(s,z) - -S"{z)p + (B'{z) - S'(z)p^ 1 (s,z))p 2 (s,z), 

z) = 2{B\z) - S'(z)p-\s, z))h(8, z), 

for q e M 2 , s > 0, and z G R, where (S,B) £ C 1 ' 1 ^). Since it is unknown whether 
■0 S C 3 (D), equation (|4.6p holds in the weak sense. It follows from (|4.4|) that 

|Ii,-(VVi,^)| < A(e), 

where A depends only on e. Furthermore, w satisfies the following boundary condi- 
tions: 

u) = on 12 = 0, 1. 
As usual for energy estimates, let r] be a C^°-function satisfying 

?? = 1 for \s\ < L, r] = for \s\ > L + 1, |t/(s)| < 2. (4.7) 
Multiply t] 2 (xi)uj and integrate it on both sides of (|4.6p . then integrate the left side, 

1 'j - 

7? 2 |Vw| 2 

p(|v^| 2 ,^)' 



and plug the explicit forms of Ay, pi(| VV>| 2 , V'); an d />2 ( I V-j/j | 2 , -0) into it. We obtain 



dxidx 2 = y^k, (4.8) 

_D n — ^ 



where 



/ 5(|V^| 2 ,^(p 2 (|V^| 2 ,^) C 2 -|VVi| 2 ) Xl X2 ' 



ll D ¥mm^ dMxid ^ 

P2(\V^\ 2 ,^)V^-Vy 2 , , 

- rjuj d Xl dx 2 , 

p 2 (\Vif)\ 2 ,tjj) 

p 2 {\V$\ 2 ,iP)c 2 - |V-0| 2 



rju} 2 dxidx<2 



h = - [[ {B"$)p{s,i>) - -S"^)p^) V Wdx 1 dX2, 

J Jd 7 



p 2 {b'^)-s'{^-\\vm 2 ^)Y „ 

d - 2 ^.ii9 ^.9 ^TjB — ^(|VV| w cb^ 



p 2 (|V</f,?/>)c 2 - |V# 
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Now we make the estimates. First, by the Holder inequality, it is easy to see that 

h + h + h< 0. 

Second, since ||(5 — S_,B — Q^c x < l ([Q,\\) — $ anci m ^ we have 

Thus, 

fL+l rl 

\h\ < OS 1 "? / / uj 2 dxidx 2 , (4.9) 



'-L-1 JO 

and p < q(D;x2), where C is independent of e. Thus, from (|4. 8[) and the definition 
of 77, if 82 is sufficiently small, we obtain 

-L r l 

\V uj\ 2 dx2dxi 

-L JO 
< \I 2 +h\ + \h\ 

<C(e)U +1 )( f {\Vuj\ 2 + \Vu\+u 2 )dx 2 )dx l 

J — L — 1 " L " 

+C , 4~' 3 f [ \Vu}\ 2 dx 2 dx 1 

J-L Jo 

<C(e)( [ L + j L+l )n'\\Vio\ 2 +Lo 2 )dx 2 )d Xl 



, L f\ ,2 

+ 7: / / dx 2 dx±. 



L JO 



Notice that 00 = on X2 = 0, 1. It follows from the Poincare inequality that there 
exists a constant C independent of I such that 

J j^\Vu\ 2 dx2dx l <c( K j +J + |Vw| 2 dx 2 )(ixi (4.10) 

for large L. It follows from (I4.5D that 

L /•! /•— L rL+l rl 



Vuj\ z dx 2 dxi < I / +/ J(/ |Vwrd2C 2 )da;i < C 

for some uniform constant C independent of L and for some constant C. Passing 
the limit L — > 00 yields 

f 00 rl 

/ I Vu\ 2 dx2dx\ < C. 
J-00 Jo 

Hence, 

/ + / )(/ l Vw | 2dx2 ) rfxi ^ asL-^oo. (4.11) 
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Using (|4.10p by passing the limit I — > oo as before again, we obtain 

j | Vu)\ 2 dx2dxi = 0, 

oo JO 

which implies uj = 0. Therefore, 

which solves the following boundary value problem: 

HO) = 0, (4.12) 

= m, 

which completes the first part. 

3. Explicit form of •tp(x2). Suppose that there are two solutions tpi and ^2 to 
(|4.12p . Let (j) = ■ipi — , ip2- Then (p satisfies 

\0(O) = 0(1) = O, 1 • j 

where 

_ = r 1 j 5(|^| 2 ,^)-2 / 5 1 (|^| 2 ,^)|^| 2 ^ - = f 1 -him 2 ,^' dc 

tffl^'l 2 ,^')^, <?= / ed^'l 2 ,^')^, 

JO 

with ?/> = st/>i + (1 — s)tp2, where •& and are defined in (|4.6p . Multiplying <^> on 
both sides of equation (|4.13j) and integrating it over [0, 1], we have 

(aft 2 + b<j>'(j>)dx2 = [ (c<f>'4> + d4> 2 )dx 2 . 
jo Jo 

Thus, we have 

-I r i Aft 

^^dsdx2 
Jo p 2 (|V' / | 2 ,'0) 

1 f 1 (B"$)p{\i>'\ 2 ^) - -S"^)f(i>' 2 ^))4> 2 dsdx 2 
Jo 7 

7o 7o 2£p - (7 + l)5p7 7 y 2£p - (7 + 1)5^7 

1 " Vfif-^f*^ 



7o 2£p 2 - (7 + 1)5/5T+ 1 



24 GUI-QIANG CHEN XUEMEI DENG WEI XIANG 

The sum of the last three terms is negative. By the smallness of 5 and the Poincare 
inequality as in Step 2, we have 

-l 

(j}'\ 2 dx 2 < 0, 







which yields = 0. Thus, the solution to (|4.12p is unique. Obviously, we know that 

^ = ip(x 2 ) = / p (s;p )u (s)ds. 



is a solution to the boundary value problem (|4,12p . In fact, from (|2.19p . 



2 

we have 



+ S(ip(x 2 ))pl 1 (x 2 ;p ) = B(tp(x 2 )), po(x 2 ;po) = ( j^T^T^y ) 



- u' (x 2 ) = B'{ifj)p (x 2 ;p ) - -S'(ip)p2(x 2 ;po) 



,Po{x 2 ;po)J X2 7 
This completes the proof. □ 

It follows from Lemma EJ] that the flow induced by the stream function ip satisfies 
(ll.20p - (jl.2ip in the upstream. The similar properties in the downstream can be 
obtained in the same way. 

As indicated at the beginning of this section, an important maximum estimate 
for the stream function can be yielded as a consequence of the far field behaviors, 
and one could see the reason for the way defining B in (13. ip . 

Proposition 4.2. Suppose that (|1.17p holds with 5 < minj^i, 5 2 }, and (S(x 2 ), B{x 2 )) 
satisfy (I1.18|) . Then the solution to Problem 3 satisfies 

0<ip <m inn. (4.14) 
Proof. It follows from Proposition 14 . 1 1 that 

f-X2 

ip(x\,x 2 ) — > / po(s;po)uo(s)ds uniformly as x\ — > — oo, 
Jo 

and 

Y>(xi, x 2 ) — > / p\(s;pi)ui(s)ds uniformly as x i — >■ oo. 
Jo 

Therefore, for any e > 0, there exists L > such that 

— e < ip(x±, x 2 ) < m + e if \xx\ > L. (4-15) 

We claim that 

— e < il){x\,x 2 ) < m + e \i\xx\<L. (4-16) 

We show our claim by contradiction. More precisely, suppose that there exists a 
point X max = {xiq,x 2 q) with |xio| < L such that 

VK^max) = max ip(xi,x 2 ) > m + e. 

XG{|a;i|<L} 
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Let p = p(ip(X max )). We have 

iy(0,V^max))<%VX^max) = f>(B' - ±S' p^ 1 ) 

= p( B '-i§s'mx max )) 

where we have used the fact that Vip(X max ) = 0, hence /5 7_1 = j(^(I max )). Thus, 
we have 

< i 4J (0,^(X max ))%V(^max) < 0, 

which is a contradiction. That is, 

— e < ip(xi, X2) < m + e in {ip > m} f]{\x\\ < L}. 
Since ln(5 _7 ,6) < in the domain {ip < 0}, we can similarly show that 

-e < ip(x 1 ,x 2 ) <m + e in {ip < 0} n {\x\\ < L}. 
Combining these estimates together, we obtain 

— e < ip(xi,x 2 ) < m + e in CI. 
Since e is arbitrary, we have 

< ip{x\,X2) <m in 0. 
This completes the proof. □ 

5. Existence and Uniqueness of Solutions of Problem 2 

Propositions 13.11 and 14.21 imply that the solutions established in Proposition 13.11 
are the solutions of the boundary value problem (|2.33[) . Problem 2. 

Proposition 5.1. Let the boundary dCl satisfy (|1.6p - (|1.9p . Let (|1.17p hold with 
8 < min{<5i, 82}, and (S (X2) , B (X2)) satisfy (|1 . 18[) . Then there exists < 5\ < 
min{<5o,5i} such that, if \\(S — S_,B — fi)||ci,i([o,i]) < o~ < <5 < <5i and m € 
, mi) with mi = < m, where m is defined in Proposition 2.2, then Problem 

2 (i.e. (|2.33p ) has a uniformly subsonic solution ip G C 2 ' a (Cl) satisfying 

< ip < rn in Cl. 

Next, we will use the energy estimates again to show that uniformly subsonic 
solutions of Problem 2 are unique. 

Proposition 5.2. Let the boundary dCl satisfy (|1.6p ~ (jl.9p . TTien i/iere exists 83 G 
(0, <5o] swc/i i/tai, «/ 

(i) \\(S-S,B-B)\\cia {[ o,i])<8 with0<5<8 3 , 

(ii) m G (0 73 , m), 

then there exists at most one solution ip of Problem 2 satisfying 

0<ip(x 1 ,x 2 ) <m, |VV>| 2 < S(e) - 2e /or some e > 0. (5.1) 
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Proof. As before, let ip\ and V>2 be two solutions to (|2.33p . Set ip = tpi — tp2 ■ Then 
ip satisfies 

{di(aijdjip) + di(biip) = Cid^ + dip in 0, 
^ = onWiUW 2 , 



(5.2) 



where 



i j 



1 AjjjDtJ) dg b . = _ f 1 P2(\V4>\ 2 ,4>)dd dc 
o p 2 (\V^\ 2 ^) S ' Jo p 2 (\V$\ 2 ,$) S ' 

0(|V^| 2 ,$)$$ds, d= f 0(\Vip\ 2 ,ip)ds, 



ip = sipi + (1 — s)ip2, A-iji © and t? are defined as (|4.6p . except we replace (S,B,p) 
by(S,B,p). 

Multiplying rj 2 ip + and integrating on both sides of (|5.2p . where 77 is defined in 
(|4.7p and ■0 + = max{ , 0(x), 0}, then, similar to the proof of Lemma 14.11 we have 

|V^| 2 dxida; 2 < C(B,e) [ [ |V^| 2 dxidx 2 . 
nn{\xi\<i}n{4>>0} J J nn{i<\xi\<i+i}n{i>>0} 

Since the solutions ip\ and ip2 have the same far field behavior, and note that 
and |V?/>| — > as \x\\ — ¥ oo, we have 



// 



\Vipfdx 1 dx 2 = 0, 

nn{i>>o} 
Similarly, we can show that 

\Vip\ 2 dx 1 dx 2 = 0, 



// 

■J Jnt 



>nn{ip<o} 

which implies that ip = 0. This completes the proof. □ 

6. Refined Properties of Stream Functions for Problem 1 

In this section, we derive some refined properties for solutions to the boundary 
value problem (|2.33p . Problem 2. More precisely, it is shown that ip X2 is always 
positive, together with the asymptotic behavior and the estimates obtained in §3- 
5, yields that (p,u,v,p) induced by ip satisfies the original Euler equations, the 
boundary conditions, the constrains on the mass flux, the Bernoulli constant, and 
the entropy equation. 

Lemma 6.1. Let the boundary d£l satisfies (|1.6p - (jl.9p . Then there exists 84 E (0, 5o] 
such that, if 

(i) \\(S -S,B -B)\\ c i.i m]) < 5 with < 5 < 64, 

(ii) m G (5? ,fh), 

(iii) ip satisfies (|5.ip and solves Problem 2, 
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then ip satisfies 

< ip < m intt, (6.1) 

and 

ip X2 > inCl. (6.2) 

Proof. We divide the proof into four steps. 

1. Equation and the boundary condition. From (|5.ip and the boundary conditions: 
V'(O) = and = m, we have 

%1> X2 > on dU. (6.3) 

Let w = ip X2 . From Lemma 14.11 w satisfies 

Ay(V^) a(P2 {\Vljj\ 2 ^)d i i> w 



9i \p2(\vi>\ 2 ,i>) 9jW ) dl \ p 2 (|v^| 2 ,v) 

= @(\Vtp\ 2 ,tp)w + -&(\V'iJj\ 2 ,iJj)d i ipd i w (6.4) 

in the weak sense, where Aij, 0, and i? are defined in (|4.6p . except replacing (5, B, p) 
by (S,B,p). 

2. Positivity in Q: That is, 

io > in n. (6.5) 

First, the far field behavior of tp implies that tp X2 — > piUi > when (— l)* +1 xi — > 
oo,i = 0, 1, which implies that w(x\,x 2 ) > for > L with L sufficiently large. 
As before, multiplying (|6.4|) by w~ = min{w,0}, integrating it on both sides and 
noticing (]6.3p . we have 

dx\dx 2 



{ u<o} P 2 (lW| 2 ,</0 

< - ff (B"{i>) - -S"^)p^ l (\V^\ 2 ^))w 2 dx l dx 2 

JJ{w<0} 1 

< CS / / w 2 dx\dx 2 . 

J J{U<0} 

For each x\, we define an open set: 

K Xl := {x 2 : f(x\) <x 2 < f(x 2 ), w(xi,x 2 ) < 0} = UjZ^, 
where each P Xi is a connected open component of K X1 . Then for every x 2 £ , 

w(xi,x 2 ) = / 9a; 2 w(xi, s)(is. 

J min It, 
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Therefore, 

^ dx\ J 



w dx\dx2 = I dx\ / I w (xi,X2)dx±dx2 

{w<0} 



dxi2_] / (/ d X2 w(x\,s)ds) 2 dx2 

-l ; JP Jmin/j., 

z ^ 1 

z2i i I (d X2 w(xi, s)) 2 ds(maxP — mmP) 2 dx2 

-I i JP JminP 



Hence, 



which means 



xitK 



\Vw\ 



{w<0} 



dx\dx2 < C5 I I I Vw\ 2 dx\dx2 



{w<0} P J J{w<0} 

I Vw\ 2 dxidx2 < 0, 



<{w<0} 

Thus, ([63]) must hold. 

3. Strict positivity in £1: That is, 

ip X2 = w > in (6-6) 



for any weak solutions w to (16. 4|) . Denote w := e aX2 w, which is a nonnegative weak 
solution to 

dJ^fe^diw) + (^a - - tf(|VVf , ^)^)e ra2 ^ + Ge aX2 w = 0, 



where and are defined in (|4.6|) . and 

M 

p z \ p z p z I \ p z 

with also defined in (|4.6p . Choosing a > sufficiently large so that G > 0, then 

di^e^djw) + (^f a - ^ - m^^mije^d.w < 0. 

This implies that (16. 61) holds, so does inequality (I6.ip . 



4. Positivity on boundary. We now show that > at Wi U W2 in this step. 
Without loss of generality, we prove it on Wi- 

First, if (SB" 1 )' '(m) < 0, since tp = m on W2, then, for any (x\, f2{x\)) £ W2, 
there exists a small disk M C satisfying A/"nJ7 = (x?, f2(xi)) such that dln( ^ 7g ^ > 
in M, which implies 

Aij (Vip, ip)dijTp > in TV. 
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Combining this with the fact that ip < m, by the Hopf lemma, we have 

^ 2 (x?,/ 2 (x?)) >0. 
The remaining case is (SB f )'(m) = 0. It is easy to see that ip satisfies 
A ij (Vip,i/j)d ij ^-m) - p 2 \V(i}-m)\ 2 + R(if;-m) = 

with R = — — — P ^2 m P — " • By the Hopf lemma again, we have 

d X2 tp > in W 2 . 

Similarly, we can show that ijj X2 > on W\. This completes the proof. □ 

7. Proof of Theorem 1.1 Except the Critical Mass Flux 

We now prove Theorem 1.1 (Main Theorem of this paper), except the existence 
part of the critical mass flux which will be shown in Section 8 below. 

Let <5 := min{5i,<5 2 ,<53,<$4} > 0. If ||(5- S, B -5)|| c i,i([o,i]) < 5 with < 5 < 5 , 
for any m E (5@, 25^ 2 ), there exists a solution of Problem 2. It follows from Lemmas 
14.11 and 16.11 that the flow field induced by ip satisfies (|2.10j) . and hence Proposition 
13.11 guarantees the existence of Euler flows. Furthermore, Propositions \3. l\ and \5.2\ 
imply the uniqueness of Euler flows with asymptotic (|1.13p and (|1.16p , the mass flux 
condition (jl.lip . and the asymptotic behavior determined by (|1.20[) — ()1.25[) . 

This completes the proof. 

8. Existence of the Critical Mass Flux 

In §5-7, we have shown that, for the given Bernoulli function and the entropy 
function in the upstream satisfying (11 . 1T[) f)l . 18|) . there exists a Euler flow, as long 

as m E (<5^,2Jq^ 2 ). In this section, we find the critical mass flux, which can be 
obtained by following the arguments as in [21 El [16] . For self-containedness, we give 
the proof in this section. 

Proposition 8.1. Let the boundary d£l satisfy (|1.6p - (|1.9p . S(x2) and B(x2) satisfy 
the asymptotic condition (|1.13p and ()1.16|) for x 2 G [0, 1] respectively, and let (|1. 18[) 
hold. Then there exists rh < rn such that, if m £ (5^,rh), there exits a unique tp of 
Problem 2 satisfying 

< ip < m infl, (8.1) 

M(m) := sup - (7 - l)5(^)p 7+1 } < 0, (8.2) 

n 

where B(ip) = U( >^ + S (ip) Pq -1 (ip) . Furthermore, either M(m) — > as m — >■ rh or 
there does not exist a > such that (|2.33p has solutions for all m £ (rh, rh + a) and 

sup M(m) < 0. (8.3) 

m£(iti,m+ij) 
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Proof. For the given entropy function S and Bernoulli function B in the upstream 
satisfying (|1.13p and (|1.16p and any m G (5^,m), one can define po and 1*0(^2) as in 
§2. Note that po and uq depend on m by definition; thus in this section we denote 
them by po(m) and uo(i/j;m), respectively. 

Let {e n }^ =1 be a strictly decreasing sequence of positive numbers such that e\ < 
£o/4 and e Tl i 0. We introduce 

s if s < — 2e n , 



7Cn if«>-£n- 



Us) 

Then ( n is an increasing smooth function. We define 

A n (|W| 2 ,^;m) 

:= Cn(|V^| 2 - (7 - l)5(^)^ +1 (|V^U;m)) + (7 - 1)5(^)^ +1 (|V^|, m). 
Then there exist two positive constants A(n) and A(n) such that 

for any z G R, q G M 2 , and £ G M 2 , where 

Ig^, z; m) = /5 n (|g| 2 , *; m)5 ij - 2ft(\q\ 2 , z; m)&&. 
Thus, for any m G (5",m), there exists a solution tfj n (x;m) to the problem: 



Al 7 y(q,z;rn)d i j7p = J" n (V^,^;m) in 0, 

/2(a;i)-/i(a!i) ' 



(8.4) 



where 

j- n (v^,v;m) = {&- U'irv^Wf+p^m? 

Moreover, if 

|VV {n) | 2 -(7-l)5(V)p 7+1 (|V^U;^) < -2e n , (8.5) 
then £4 = 1. Similar to §3, we have 

< if) n (x; m) < m. 



By the definition of S and B, we can estimate 1$ in (|4.9p . which is independent of e n . 
Furthermore, it follows from Lemma 14.11 that the solution in (|8.4p satisfying (|8.5p 
has the far field behavior as (14. ljl . In addition, by Proposition 15. 2\ such a solution 
is unique in the class of solutions satisfying (|4,ip . 

Note that, in general, we do not know the uniqueness of solutions to the boundary 
value problem (|8.4p . Set 



Sn{fn) = {ip n (x;m) : ip n (x;m) solves problem (|8.4p }. (8-6) 

Define 

M n (m)= inf sup{|V^| 2 -( 7 -l)5(^ +1 (|V^|,V;H}, (8-7) 

i/> n G5„(m) q 
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and 

T n = {s : 5? < s, M n (m) < -4e„ if m £ (6?, s)}. 
It follow from Proposition 13. l\ Lemma 14.11 and Proposition 14,21 that 

[<,2f]cT„, 

hence T n ^ 0. We define m n = supT n . 

The sequence {m n } has the following properties: 

1. M n (m) is left continuous for m E (5@ ,m n ]. Indeed, let {m^ E (<5^,TO n )} and 
mffl t m - Since M n {rrin^) < — 4e n , we have 

\\^ n \x;mM)\\ c ., a <C(n). 

Therefore, there exists a subsequence i/j^(x;mh ) such that 

V; (n) (x;m^) 

Moreover, ip solves (|8.4p . and M n (m) < lim?// n )(x; mffl). Thus, 

M n (m) < -4e n . 

Note that all these solutions satisfy the far field behavior as (|4.1|) . by uniqueness of 
solutions in this class, 

M n (m) = lim^ n) (x;m^). 

2. m n < fh: If this were not true, by the definition of m n , fh G T n . It follows 
from the left-continuity of M n {m) that 

M n (fh) < -4e n . 

Thus, by means of the proof of Lemma 14, 1\ tp n (x;fh) has far field behavior as in 
(|4,ip . However, it follows from the definition of fh that 

sup{|VV (n) | 2 - (7 - l)S^)~pl +1 (\V^ n \^-m)} 

\ \po{s;fh)u (s;fh)\ 2 - (7 - l)S(s)pl +1 (s;fh), \ 

> sup max < > 
se[o,i] { \pi{s]fh) Ul (s]fh)\ 2 - (7 - l)S(s)pJ + (s;m) J 

f 2B( S )pl( S ;m)-( 1 + l)S( S )pl +1 ( S ;fn), \ 

> sup max < > 
se[0,l] { 2B{s)p 2 {s;fh)-(j + l)S{s)pJ + \s;fh) J 

> sup g\D; S )S^s)(2D(x 2 )-( 1+ l){^ly-^ 
»e[o,i] v 7 - 1 

> sup 2g 2 (D;s)S^{s)(D(x 2 ) - D) 

se[o,i] 

= 0, 

where pi{s;m) = pi(y(s);fh), u\{s;fh) = ui(y(s);fh), and y(s) is the function de- 
fined in (|2.24p . Thus, M n {fh) > 0. This is a contradiction. Therefore, m n < m. 
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Finally, {m n } is an increasing sequence, which follows from the definition of {m n } 
directly. Define 

rh = lim m r . 



n— yoo 



Then rh is well-defined and rh < m. Note that, for any m € (S^,rh), there exists 
m n > m such that M n (m) < — 4e n . Thus, 

cf) = ^ n \ X ; m) 

solves (|2.33p and 

sup { | VVf - (7 - l)S(t/0^ +1 (|Vt/> (n) I , V (n) ; m)} = M n (m) < -Ae n . 
n 

If sup M n {m) < 0, then there exists n such that 

sup M n (m) < -Ae n . 

md{&P ,m) 

Then the same argument as the proof for the left continuity of M n {m) on (8^,m n ] 
yields 

M n {rh) < -4e n . 

Suppose that there exists a > such that (|2,33p always has a solution ip for m G 

(rh, rh + a), and 

sup M(m)<0. (8.8) 

mG(m,m+cr) 

Then there exists k > such that 

sup M(m) = sup {|V^| 2 - (7 - l)5W^ +1 (|V^U;m)} < -4e n+fc . 

mG(m,m+cr) m£(m,-m+cr) 

This yields that m n+ k > rh + a, which is a contradiction. Thus, either M(m) — > 
or there does not exist a > such that (|2.33|) has a solution for all m £ (m, rh + a) 
and (18 .3p holds. This completes the proof. □ 
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